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We consider the quantum ferromagnetic transition at zero temperature in clean itinerant electron 
systems. We find that the Landau-Ginzburg- Wilson order parameter field theory breaks down since 
the electron-electron interaction leads to singular coupling constants in the Landau-Ginzburg- Wilson 
functional. These couplings generate an effective long-range interaction between the spin or order 
parameter fiuctuations of the form l/r^''"^, with d the spatial dimension. This leads to unusual 
scaling behavior at the quantum critical point in 1 < d < 3, which we determine exactly. We also 
discuss the quantum-to-classical crossover at small but finite temperatures, which is characterized 
by the appearance of multiple temperature scales. A comparison with recent results on disordered 
itinerant ferromagnets is given. 



I. INTRODUCTION 

The physics of quantum phase transitions has been a 
subject of great interest lately. In contrast to the ususal 
classical or thermal phase transitions, quantum phase 
transitions occur at zero temperature as a function of 
some nonthermal control parameter, and the fluctuations 
that drive the transition are of quantum nature rather 
than thermal in origin. Among the transitions that 
have been investigated are various metal-insulator tran- 
sitions, the superconductor-insulator transition in thin 
metal films, and a variety of magnetic phase transitions. 
Early work in this fielcB established that if the quantum 
phase transition has a classical analog at finite temper- 
ature, then the former tends to have a simpler critical 
behavior in the physical dimensions c? = 3 or d = 2 than 
the latter. The reason for this tendency is that the cou- 
pling between statics and dynamics that is inherent to 
quantum statistics problems effectively increases the di- 
mensionality of the system from d to d + z, with z the 
dynamical critical exponent. This reduces the upper crit- 
ical dimension d^, which is the dimension above which 
mean-field theory yields the exact critical behavior, by z 
from its value for the classical transition. 

One of the most obvious examples of a quantum phase 
transition, and the first one studied in detail, is the ferro- 
magnetic transition of itinerant electrons that occurs as 
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a function of the exchange coupling between the electron 
spins. In a pioneering paper, Hertaj derived a Landau- 
Ginzburg- Wilson (LGW) functional for this transition by 
considering a simple model of itinerant electrons that in- 
teract only via the exchange interaction in the particle- 
hole spin-triplet channel. Hertz analyzed this LGW func- 
tional by means of renormalization group (RG) meth- 
ods that generalize Wilson's treatment of classical phase 
transitions. He concluded that the critical behavior in the 
physical dimensions d = 3 and d = 2 is mean-field like, 
since the dynamical critical exponent z = 3 decreases 
the upper critical dimension from d+ = 4 for the clas- 
sical case to d^ = 1 in the quantum case. In order to 
study nontrivial quantum critical behavior. Hertz then 
considered a model with a magnetization that is defined 
on a space of arbitrary dimension d, while the correla- 
tion functions that determine the coefficients in the LGW 
functional are taken to be those of a 3-d Fermi gas. For 
this model, he calculated non mean-field like critical be- 
havior in d < 1 by means of a 1 — e expansion. De- 
spite the somewhat artificial nature of this model, there 
is a general belief that the qualitative features of Hertz's 
analysis, in particular the fact that there is mean-field 
like critical behavior for all d > 1, apply to real itinerant 
quantum ferromagnets as well. 

In this paper we reexamine the ferromagnetic quantum 
phase transition of itinerant electrons and show that the 
above belief is qualitatively mistaken. We first consider 
a model that is more realistic than Hertz's, viz. with an 
electron-electron interaction that is not restricted to the 
particle-hole spin-triplet channel. We find that the LGW 
approach breaks down due to the presence of soft modes 
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in addition to the critical modes, namely particle- hole ex- 
citations that couple to the conserved order parameter. 
These soft modes are integrated out in the derivation 
of the LGW functional, and this leads to singular ver- 
tices in the order parameter field theory. This is a rather 
general observation; analogous effects are expected for a 
large class of quantum phase transitions, and in general 
they invalidate the application to quantum phase tran- 
sitions of the classical LGW philosophy of deriving an 
effective logal field theory entirely in terms of the order 
parameter.l3 However, for the present problem we find 
that the mathematical difficulties generated by the sin- 
gular couplings can be handled. The resulting nonlocal 
field theory contains an effective long-range interaction 
between the order parameter fluctuations, and we are 
able to determine the critical behavior exactly for all 
d > 1. For 1 < d < 3 the result is different from ei- 
ther mean-field critical behavior or classical Heisenberg 
critical behavior. We then show that the same is true for 
Hertz's original model. Our results invalidate the 1 — e ex- 
pansion in that paper, and it corroborates jand explains 
an observation made recently by Sachdevja who noted 
that Hertz's results in cZ < 1 cannot be correct since they 
violate an exact exponent equality. 

The outline of this paper is as follows. In Sec. || we 
first define our model of itinerant electrons, and then dis- 
cuss how to derive an order parameter description for the 
ferromagnetic phase transition starting from a fermionic 
field theory. We discuss the behavior of the coefficients in 
the resulting LGW functional and show that they have 
long-range properties. In Sec. HI the quantum critical 



behavior of the resulting nonlocal field theory at zero 
temperature is determined exactly for dimensions d > 1. 
For d > 3, mean-field exponents are obtained, while for 
1 < d < 3, d-dependent exponents are found. In Sec. 
[V we investigate the behavior at small but finite tem- 
peratures which is characterized by multiple independent 
temperature scales. In Sec. ^ we reexamine Hertz's origi- 
nal model of the ferromagnetic quantum phase transition. 
We show that upon renormalization it acquires the same 
features as the more realistic model studied in Sees. || - 
[V. In Sec. VI we discuss our results, and in particular 



the relation of the present paper tojrecent work on dis- 
ordered itinerant quantum magnetsu A technical detail 
concerning perturbation theory is relegated to Appendix 
A, while in Appendix B we discuss logarithmic correc- 
tions to scaling that exist in d = 3 as well as for all 
dimensions 1 < d < 3. A short accoupt of some of our 
results has previously been published.El 



II. ORDER PARAMETER FIELD THEORY FOR 
THE QUANTUM FERROMAGNETIC 
TRANSITION 

In the first part of this section we define our model 
of itinerant interacting electrons in terms of a fermionic 



field theory. Since we are interested in a phase transi- 
tion to a magnetically ordered phase, we choose the spin 
density fluctuation as our order parameter, and proceed 
to derive an LGW or order parameter description of this 
transition. We then derive and discuss the cocflficicnts 
in this LGW functional. As we will see, the interac- 
tions in the effective LGW theory are long-ranged due to 
soft modes that have been integrated out in deriving the 
LGW functional. 



A. The Model 

The partition function of any fermionic system can be 
written in the form,El 



J D^pDxl; exp(5 [V-,?/-]) 



(2.1a) 



where ip and tp ^-re Grassmannian (i.e., anticommuting) 
fields. Dij: Dip denotes the Grassmannian functional in- 
tegration measure, and S is the action. 



S = drjdx V5'(x, r) ^ 7/.,(x, r) - dr H{t) 



(2.1b) 



Here x denotes position and t imaginary time, H{t) 
is the Hamiltonian in imaginary time representation, 
f3 — 1/T is the inverse temperature, j = 1,2 denotes 
spin labels, and summation over repeated covariant and 
contravariant spin indices is implied. Throughout this 
paper we use units such that ks = = e^ = 1. We 
start with a microscopic model of itinerant, interacting 
electrons, 



H{t) = / dx V*(x,r) 



1 

2m 



V'»(x,t) 



+ 2] '^^'^^"(^^y) '0'(x,T)'0^(y,T)'0j(y,T) 

xi^,(x,t) . (2.2) 

Here m is the electron mass, fi is the chemical potential, 
and u(x— y) is the electron-electron interaction potential. 
More realistic models to describe itinerant electron mag- 
netism including, e.g., band structure, can be considered 
along the same lines. The salient points of the present 
paper, however, are due to long-wavelength effects and 
therefore are independent of microscopic details like the 
band structure. For our purposes it therefore is sufficient 
to study the model defined in Eq. (2.2). 

In order to describe magnetism, it is convenient and 
standard practice to break the interaction part of the ac- 
tion S, which we denote by Sint, into particle- hole spin- 
singlet, particle-hole spin-triplet and particle-particle or 
Cooper channel contributions, which we denote by S^^j, 
respectively.0 For simplicity, we assume 



^^^^ and <?'^'') 
^int "^II^ '-'mt ; 
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that the interactions are short-ranged in all of these chan- 
nels. In a metallic system this is justified due to screen- 
ing, and an effective model with a short-ranged inter- 
actions carube derived starting from a bare Coulomb 
interaction.Q The spin-triplet interaction warrants spe- 
cial attention, since the interactions between spin density 
fluctuations that are described by S*/^] are what causes 
ferromagnetism. We therefore consider this part of the 
action separately, and write 



S = So + S. 



it) 



with 



,(t) 
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J dx dr ns (x, t) • (x, t) 



(2.3) 



(2.4a) 



where is the electron spin density vector with compo- 
nents, 



<(x,T) = i V5^(x,T)a,'^.V'^(x,T) 



(2.4b) 



Here the ct" (a = 1,2,3) are the Pauli matrices, and Ft 
is the spin-triplet interaction amplitude. 5*0 in Eq. ( |2.3| ) 
contains the free electron part and all interaction parts 
other than the particle-hole spin-triplet contribution to 
the action. It reads explicitly. 



Sn = 



dr / dx 



d 

OT 
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~^ 
rP 

-Tc / dr 
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rix,T) [-VV2m-Ai] V,(x,t) 
J '^'^ J c^x n(x, t) n(x, r) 
dx nl^{x, t) nc,ij{x,T) 



(2.5a) 



Here n is the electron number or charge density, 

n(x,r) = V;''(x,T)V^,(x,r) , (2.5b) 
and Tic.ij is the anomalous or Cooper channel density, 



(x,t) = {I - 6,j) ■)P,{x,t)^j{x,t) 



(2.5c) 



(no summation). To avoid double counting, only the 
long-wavelength fluctuations of the charge density and 
t he an omalous density must be taken into account in Eq. 
(2.5a). a Fs and Fc are the interaction amplitudes in the 
particle-hole spin-singlet and in the Cooper channel, re- 
spectively. For our purposes it is sufficient to treat Ft, 
Fs, and Fc simply as numbers, although in a more com- 
plete model thej5£^are complicated short-ranged interac- 
tion potentialsEfl For later reference we note that they 
result from averaging the original interaction over differ- 
ent regions of phase space, so that the three interaction 
amplitudes are independent parameters. 



B. Order parameter field theory 

The standard theory of continuous thermal phase tran- 
sitions proceeds from a microscopic model by identify- 
ing the order parameter relevant for the transition, and 
by integrating out all degrees of freedom except for the 
long wavelength order parameter fluctuations, or criti- 
cal modes. The result of this procedure is a LGW the^. 
ory, i.e. an effective field theory for the critical modes.EHl 
The rationale behind this approach is that these fluctu- 
ations, which are slowly varying in space, determine the 
behavior near the critical point. The same basic idea 
has been applied to quantum phase transitions, with the 
only principal difference being that the critical modes 
are now slowly varying in both space and (imaginary) 
time. While we will use this approach here, motivated in. 
part by previous work on itinerant electronic systems,u 
we mention that in general one should worry about both 
the critical modes, and all other slow or soft modes, even 
if these other soft modes are not 'critical' in the sense 
that they change their character at the phase transition. 
Of course this concern is not restricted to quantum phase 
transitions; additional soft modes can and do occur also 
at thermal phase transitions. However, while this is an 
exceptional occurence in the case of the latter, we will 
argue below that it poses a more serious problem for the 
former, since at T = there are more soft modes than at 
finite temperature. We will see that in the present prob- 
lem, and in a large class of other quantum phase transi- 
tions, such additional modes are indeed present and lead 
to complications if one insists on a description entirely 
in terms of the order parameter. For the problem under 
consideration, however, we will be able to overcome these 
problems and determine the critical behavior within the 
framework of an LGW theory. 

The techniques for deriving an order para mete r field 
theory, starting with Eqs. (2J) through (|2.5|), are 
standard.El We introduce a classical vector field M(x, r) 
with components M°, the average of one of which {M^, 
say) is proportional to the magnetization m, and de- 
couple the four fermion term in 5*^^^] by performing a 
Hubbard-Stratonovich transformation. Subsequently, all 
degrees of freedom other than M are integrated out. This 
procedure in particular integrates out soft particle-hole 
excitations, which are the additional soft modes men- 
tioned above. We obtain the partition function Z in the 
form. 



^ = e--Fo/T / ^j]y[j cxp(-$[M]) 



(2.6a) 



with Fq the noncritical part of the free energy. The LGW 
functional <I> reads, 



$[M] = Y J '^^ ■ 



In ( exp 



-Tt J dx M{x) ■ ns{x) 



(2.6b) 



So 



3 



Here we use a four- vector notation with x = (x, t 
r/3 



and 



.0 



J dx = J dxjg dr, and {■ ■ ■)so denotes an average taken 
with the action Sq. A formal expansion of <f> in powers 
of M takes the form, 



$[M] = i / dxdy 



S{x - y) - X^'^\x - y) 



xM{x) ■ M{y) 



n—3 



dxi 



dXr, 



XXai?..a„ ("^1 2;„ , . . . , a;„_i Xn) 

xM'^'ixi) ... M'^-ixn) , 



(2.7a) 



where 6„ — (— l)"/n! and w e ha ve scaled M with T^^. 
The coefficients x*^"-' in Eq. (2.7a) are connected n-point 
spin density correlation functions of a reference system 
whose action is given by Sq. They are defined as, 

xi:LM,---,^n) = {n:H^i)---<"ixn))'s, , (2.7b) 

where the superscript 'c' denotes a cumulant or con- 
nected correlation function. The tensor structure of 
X*-"' is restricted by rotational invariance in spin space. 
For instance, x^'^'' is proportional to the completely 
antisymmetric third-rank tensor e"^'"^. For our simple 
model, the reference ensemble Sq consists of free electrons 
with short-ranged particle-hole spin-singlet and particle- 
particle model interactions. As mentioned in Sec. II A 
above, the model can be made more realistic, if desired, 
by e.g. including band structure effects. To this end 
one would simply replace the x*-"-* above with correlation 
functions for band electrons. 



C. The coefficients of the LGW functional 

To proceed, we have to calculate the coefficients in 
the LGW functional, Eq. ( 2.7a ), i.e. the connected spin 
density correlation functions of the reference system 5*0. 
Normally, one would localize the individual terms in Eq. 
( p. 7a ) about a single point in space and time, and expand 
the correlation functions in powers of gradients and fre- 
quencies. However, for the system under consideration 
such a gradient expansion does not exist. To illustrate 
this point, let us first consider the spin susceptibility 

(2) 

x\ij = SabXsi whose Fourier transform is defined by 

X.(q,aj = J rf(xi -X2)d(Ti -ra) e''^<^'-^-'> 

xe'^-^^'-^-^ Xs{xi - X2) . (2.8) 

where fin = ^irTn is a bosonic Matsubara frequency. In 
order to study the critical behavior we have to determine 
X^^^ in the long- wavelength and low- frequency limit. In a 
system with a conserved order parameter, the frequency 
must be taken to zero before the wavenumber, otherwise 



one never reaches criticalityJlil In the critical region we 
therefore have |f2„| << q. In this limit the spin suscep- 
tibility has the structure 



Xs{q,^n) = xo{q) [i-|aj/9 + 



(2.9) 



where q and f2„ are being measured in suitable units, 
and xo (?) is the static spin susceptibility of the reference 
system. 

The static spin susceptibility of non-interacting elec- 
trons, i.e. the Lindhard function, as a function of the 
wavenumber q is analytic at g = 0. However, for any 
electron system with a nonvanishing particle-hole spin- 
triplet interaction. Ft 0, there is a-ponanalytic correc- 
tion to the static spin susceptibilityJla Although our bare 
reference system 5*0 does not contain such an interaction 
amplitude, in a generic reference system a nonvanishing 
Ft is generated in perturbation theory. For being generic 
in this sense, it suffices that the bare system has a non- 
vanishing Cooper channel interaction amplitude Fc 0. 
This is shown in Appendix A. Effectively, we therefore 
have to calculate the static spin susceptibility of a para- 
magnetic Fermi liquid with a nonvanishing Fj. It has 
been shown in Ref. |lj that the result is, 

Xo(g^O) = Co- Cd^K/'^-i-ca + . , (2.10a) 

whB£e C2 > 0, and Cd-i is proportional to F^ for small 
Ft.t3 The nontrivial, and for ou r purposes most interest- 
ing, contribution in Eq. ( p. 10a ) is the nonanalytic term 
^ q'^~^. The above form holds for 1 < d < 3. In d = 3 
the nonanalyticity is of the form In q, and for d > 3 
the leading (/-dependence is given by the quadratic term. 

This feature of a Fermi liquid is of crucial importance 
for what follows. It therefore warrants some discus- 
sion, fijien though it has recently been discussed in great 
detail.li3 It is well known that in a Fermi liquid, both in 
d — 3 and in d ~ 2, the specific heat and the quasipartide 
lifetime are nonanalytic functions of the temperature.tj 
The physical reason for these effects is, in the language of 
the present paper, the soft particle-hole excitations that 
are always present in a electron system at zero tempera- 
ture. Although they become massive at nonzero temper- 
ature (see Eq. ( p.lOb| ) below), this is sufficient to make 
various observables nonanalytic functions of temperature 
at T = 0. The nonanalytic wavenumber dependence of 
the spin susceptibility at T = is just another mani- 
festation of this effect. It is analogous to a feature of 
disordered electron systems,rJt)r which it is known that 



Xoiq ^ 0) ~ const. -I- g'' 



Since the physical reason 



for this effect is again the presence of soft modes, it is 
not qualitatively tied to disordered systems. The only 
difference between disordered and clean systems in this 
respect is the nature of the soft modes: In the former 
they are diffusive, leading to an exponent d — 2, while in 
the latter they are ballistic and lead to an exponent d — 1. 

In the present context, the nonanalytic behavior of 
Xo implies that the standard gradient expansion men- 
tioned at the beginning of this subsection does not ex- 
ist. Rather, the interaction between the order parameter 
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fluctuations is effectively of long range, and in real space 
takes the form r~'-^'^~^\ This is true only at T = 0. At 
finite temperature, where one has to perform a frequency 
sum rather than a frequency integral to calculate the cor- 
relation function, the nonanalytic. term is replaced by a 
term of the schematic structureJlj 



(d-i) 



(2.10b) 



(2.11b) 
and ( ^.llap 



( |2.7a| ), contains a 



so for fixed T > an analytic expansion about q — ex- 
ists, and the standard local LGW functional is obtained. 

We now turn to the higher spin density correlation 
functions x''"''- The same physics that causes the non- 
analyticity in the spin susceptibility, Eq. ( p. 10a ), leads 
to an even stronger effect in the higher correlation func- 
tions, and results in x^"^ for d < n — 1 not being finite 
in the limit of zero frequencies and wave numbers. One 
finds that in this limit x*^") is schematically given by 

X^") ~ w„ + , (2.11a) 

where u„ and w„ are finite numbers and we have cut off 
the infrared divergence by means of a cut-off momentum 
p. Rotational symmetry in spin space requires U2n+i = 
0. Again, the nonanalytic behavior is confined to zero 
temperature, at finite temperature x''"'' is given by 

which is finite for p 0. Equations 
imply that our LGW functional, Eq, 
nonanalyticity which has the form of a power series in 
M/p. In order to specify the LGW functional, we still 
need a physical interpretation of the infrared cutoff mo- 
mentum p. This will be given in the next section. 

To conclude this section we show that the static cor- 
relation functions discussed above provide more impor- 
tant nonanalyticities than higher order terms in a fre- 
quency expansion. To see this, we anticipate a result 
from the next section. We will discuss a fixed point of 
the renornialization group transformation where the dy- 
namical exponent is z = d, i.e. rj„ scales like q'^. Now let 
us look at an expansion of x^"^ in powers of frequency. 
In our clean system with ballistic modes the term of m*'' 
order in f2„ carries an additional factor of in the de- 
nominator, as can be seen in Eq. (2.9) for the spin sus- 
ceptibility. Thus, if the term of zeroth order in the fre- 
quency scales like q'^+^-^ near criticality, then the mth 
order term scales like qd+i-n+m(d-i) _ Yqi d > 1 the 
static susceptibility therefore has a stronger divergence 
than the frequency corrections. However, in d = 1 the 
leading divergence is provided by a frequency dependent 
term, as will be discussed in Sec. below. 



theory, and by studying the renormalization group prop- 
erties of the Gaussian fixed point. Wc then analyze the 
non-Gaussian terms in the field theory, and show that 
they are irrelevant, in the renormalization group sense, 
with respect to the Gaussian fixed point for all dimen- 
sions d > 1. This implies that the Gaussian theory yields 
the exact critical behavior for all of these dimensions, ex- 
cept for logarithmic corrections to scaling in d = 3 that 
are discussed in Appendix B. We then construct the 
equation of state near the critical point. This requires a 
more detailed discussion of the irrelevant non-Gaussian 
terms in the field theory, since the equation of state is 
determined in part by dangerous irrelevant variables. 



A. The Gaussian fixed point 



According to Eqs. (g/7|), (^, and ( |2.10a| ), the Gaus- 
sian part of the LGW functional <I>[M] is, 

$2[M]=i^^M(q,a„) [io + fld-ig'^-^ 
q a„ 

+ a2q^+ an \n„\/q] ■ M(-q, -i},,) , (3.1a) 



where 



to = l~Tt Xs(q^O,n„ = 0) 



(3.1b) 



is the bare distance from the critical point, and a^-i, 0.2, 
and ao are positive constants. 

W e first analyze the critical behavior implied by Eqs. 
( |3.l[ ). Later we will show that for d > 1 fluctuations 
are irrelevant, and the critical behavior found this way 
is exact for these dimensions . Fo ur critical exponents 
can be directly read off Eq. ( 3.1a ). These are the cor- 
relation length exponent v, defined by ^ ~ t~^, with ^ 
the correlation length and t the (rcnormalized) dimen- 
sionless distance from the critical point; the exponent 77 
that determines the wavenumber dependence of the mag- 
netic (i.e., order parameter) susceptibility at criticality, 
(Afa(q,0) Ma(— q, 0)) ~ q~^~^^; the dynamical scahng ex- 
ponent z that characterizes critical slowing down by re- 
lating the divergence of the relaxation time, Tr, to that of 
the correlation length, ^ and the exponent 7 that 
describes the i-dependence of the static magnetic suscep- 
tibility Xrn_= (Ma(0,0)Ma(0,0)) i-T. An inspection 



of Eq. (3.1a) yields. 



f l/(d- 1) for 1 < d < 3 
\l/2 ford>3 



(3.2a) 



III. THE CRITICAL BEHAVIOR AT ZERO 
TEMPERATURE 



3-d 
.0 



for 1 < d < 3 
for d > 3 



(3.2b) 



Here we discuss the quantum critical behavior of the 
Gaussian part of the LGW theory defined in the last sec- 
tion. We do so both by exphcitly solving the Gaussian 



3 for d > 3 



(3.2c) 
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7=1 for all d > 1 



(3. 2d) 



In d = 3 there are logarithmic corrections to scaling, see 
Appendix B. 

For later reference, we also discuss the critical behav- 
ior given by Eqs. 



(3.1 



L) and (3.2) from a renormalization 
group point of view. Let b be the renormalization group 
length rescaling factor. Under renormalization, all quan- 
tities change according to A ^ A{b) — b^-^^ A, with [A] 
the scale dimension of A. The scale dimension of the 
order parameter is. 



[M(q,aO] 



-l + ?7/2 



or, equivalently. 



[M(x,r)] = (d+l)/2 



(3.3a) 



(3.3b) 



At the critical fixed point, an and either a^-i (for 1 < 
d < 3), or 02 (for d > 3) are not renormalized, i.e. their 
scale dimensions are zero. Using this, and [q] — 1, \il 



z immediately yields 77 and z as given by Eqs. ( ^ , |3!2c| ). 
Equation (3.2a) follows from the relevance of to, or its 
renormalized counterpart, t, at the critical fixed point. 
That is, the scale dimension of t is positive and given by 

=[t]=2-TJ. 



B. The non-Gaussian terms 

We now show that all of the non-Gaussian terms in the 
field theory are renormalization group irrelevant with re- 
spect to the Gaussian fixed point discussed in the last 
subsection. To do this we determine the scale dimen- 
sions of the coefficients m„ and u„ of the higher-order 
term s in th e LGW functional which had been defined in 
Eq. (2.11a). Since the latter have been defined in Fourier 



space, we first take the Fourier transform of the n sum- 



mand in Eq. (2.7a). This yields 



--(n-2)[T/t/]-n[M(q,aO] 
-(n-2){d + z) + n{l-r]/2) , 



which is always negative. For odd n, and for d < n — 1 
for even n, we need to consider, 



-(n- 2)(d- l)/2 
2-n((i-l)/2 



for 1 < d < 3 
for d < 3 



(3.5b) 



(3.4) 



We see that all of the u„ are irrelevant for d > 1, while 
in d = 1 they all become marginal. We conclude that for 
d > 1 the critical Gaussian fixed point is stable, and so 
the exponents given in Eqs. (3.2) are exact. In d = 1 an 
infinite number of operators seems to become marginal, 
so naively one would conclude that the upper critical di- 
mension is given by d^ = 1. However, one has to keep in 
mind that the functional form of correlation functions in 
d = 1 can be qualitatively different from that in d > 1, so 
our power counting may be valid only for d > 1, and the 
scale dimension of some operators may change discontin- 
uously from irrelevant in d > 1 to relevant in d = 1 . This 
actually happens, as will be discussed in Sec. ^ below. 



C. The equation of state 

We now determine the equation of state, and calcu- 
late the critical behavior of the magnetization and the 
magnetic susceptibility. Since we have shown in the last 
subsection that fluctuations are irrelevant for the critical 
behavior for all d > 1, we can determine the equation of 
state by simply calculating the saddle point contribution 
to the free energy. However, in order to do so, we have 
to include the higher order terms in the LGW functional. 
Although they are irrelevant operators for d > 1, they are 
potentially dangerously irrelevantEs with respect to the 
magnetization m. We will see below that all of the 
are indeed dangerouly irrelevant, since to is a singular 
function of the Vn for w„ ^ 0. 

We determine the mean field equation of state by cal- 
culating the saddle point contribution to the free energy. 
To do so, we replace the o rder parameter field M(a;) in 
the LGW functional, Eq. ( ^.7a ), by the average magne- 
tization m and determine the stationary point of $ with 
respect to to. This yields the equation of state in the 
form 



where we have used Eq. (3.3a). Next we need to assign a 
scale dimension to the cutoff momentum p in Eq. (2.11a). 
The most obvious guess from a scaling point of view is to 
identify p with the inverse correlation length which 
makes [p] — 1. We will ascertain in Sec. IIIC below that 



t his is indeed the correct choice. It then follows from Eq. 



( 2.11a ) that of the two parameters u„ and w„ , Vn has the 



larger scale dimension, and hence is more relevant, for 
d < n — 1, while u„ is more relevant for d > n — 1. For 
even n (remember that U2n+i = 0) and d > n — 1 the 
most relevant parameter thus has a scale dimension 



\U2n 



-(d+l)(n-l) + 2 



(3.5a) 



tin + U4 + TO** Vn {m/pY 



= H 



(3.6) 



n— 3 



where an external magnetic field H was added. We have 
only kept the leading terms in each order of to, and we 
have suppressed all numerical prefactors since they are 
unimportant for our purposes. The equation of states 
now has the form of a power series in m/p. This implies 
that the cut-off momentum p effectively scales like to. 
Thus, all higher order terms effectively have the same 
power of to, viz. m'^, and the equation of states now 
reads 



tm + U4 to +1)171 = H 



(3.7) 
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with U4 fro m E g. ( ^.lla ) and v another finite coefficient. 
From Eq. (3/7) we immediately obtain the exponents /3 
and 5, defined by m{t, H = 0) ^ tf^ , m{t = 0, i/) - H^^\ 



as 



forl<d<3 
1/2 ford>3 

^ (d ior l<d<3 
I 3 for d > 3 



(3.8a) 



(3.8b) 



In d = 3 logarithmic corrections to scaling occur, see 
Appendix B. 

We are now in a position to determine the exact scale 
dimension of the cut-off momentum p and to verify the 
idetification of p with made in the last subsection. 
As we have seen after Eq. (3.6), the cutoff p scales like 
m, so that [p] = [m]eff, where [rnjcff = P/i^ is the effective 
scale dimension of m, i.e. the scale dimension after the 
effects of the dangerous irrelev ant v a riable s have been 
taken into account. From Eqs. ( 3.8a ), ( 3.2a| ) we see that 
f3 ~ v, hence [m]off = 1. Consequentely, [p] = 1 which 
justifies the identification p with for scaling purposes 
made in Sec. IIIB. These results complete the proof of 



the statement that the system is above its upper critical 
dimension for d > I. 



IV. BEHAVIOR AT FINITE TEMPERATURES, 
AND THE QUANTUM TO CLASSICAL 
CROSSOVER 

In the first part of this section we discuss various 
sources of temperature dependence in our field theory 
and the corresponding scaling behavior of the temper- 
ature. In the following two subsections we apply the 
results to a calculation of the equation of state at finite 
temperatures and of the specific heat, and we discuss 
the crossover from the quantum Gaussian fixed point to 
the classical Gaussian fixed point. In the last subsec- 
tion we discuss the analogous temperature dependence 
in systems with quenched disorder. 



A. Scaling behavior of the temperature 

T he tem perature dependence of the LGW functional, 
Eq. ( 2.7a| ), is due to a number of entirely different effects. 
First, the spin-density correlation functions x'"'' that de- 
termine the coeffients of the LGW functional are tem- 



perature dependent as given in Eqs. (2.10b) and (2.11b) 



These are correlation functions for the reference ensem- 
ble that is far from any critical point. Therefore, the fre- 
quency or temperature in x'") scales like the momentum, 
and consequentely we have a temperature scale whose 
scale dimension is, 



This we will refer to as the ballistic temperature scale. 

A second source of temperature dependence is the 
usual dynamical scaling. This originates from the fact 
t hat the time integration in the LGW functional, Eq. 
( ^.7a ), extends over the finite interval [0, 1/T]. Dynam- 



ical scaling is hence equivalent to finite-size scaling in 
the direction of imaginary time. The scaling behavior 
of the temperature due to dynamical scaling is described 
by the dynamical exponent z as given in Eq. (B.2c). We 
therefore have a second temperature scale whose scale 
dimension is 



[rjcrit — 



(4.1b) 



which below we will refer to as the critical temperature 
scale. 

Finally, the distance from the critical point, t, is tem- 
perature dependent. For the bare to, this is simply the 
usual dependence that is familiar from Fermi liquid 
theory. The scale dimension of this Fermi liquid temper- 
ature scale is thus 



[T] 



FL 



l/2iy 



(4.1c) 



However, upon renormalization t acquires a more intri- 
cate temperature dependence. Since the loop or fiuctu- 
ation corrections involve integrals over critical propaga- 
tors, the latter depends on the critical temperature scale. 
However, it also depend s on t he scale dimensions of the 
vertices m„ and w„, Eq. (2.11a), that are dangerous irrel- 
evant variables with respect to the T-dependence of tha 
magnetizatiOjtL.and the magnetic susceptibility. Millistll 
and SachdevES have shown that the resulting tempera- 
ture scale has a scale dimension 



[Thud = z/{l ~ vO) 



(4.1d) 



[r]baU = 1 



(4.1a) 



where 6 is the scale dimension of the appropriate danger- 
ous irrelevant operator, U4 in our case. This we call the 
fiuctuation temperature scale. 

Of all the temperature scales that are present in a given 
quantity, the one with the largest scale dimension will be 
the dominant one. This means, for instance, that the 
Fermi liquid scale will be subdominant compared to the 
ballistic one for all 1 < d < 3; for d > 3 the two scales 
are indistinguishable. Since z > I for all d > 1 it also 
means that the critical temperature scale will in general 
be dominant over both the ballistic scale and the fiuctua- 
tion scale. There are, however, two possible mechanisms 
that can invalidate this conclusion. The first possibil- 
ity is that some quantities do not depend on the critical 
temperature scale. For those the leading temperature 
dependence will be given by cither the ballistic tempera- 
ture scale, or by the fiuctuation scale, depending on the 
values of [Tjbaii and [Tjauct- As we will see, this possi- 
bility is realized for the magnetization and the magnetic 
susceptibility. The second possibility is that a subdom- 
inant temperature scale is dangerously irrelevant. This 
possibility cannot be ruled out by general scaling consid- 
erations. 
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B. Magnetization, and magnetic susceptibility 

Since we work above an upper critical dimensional- 
ity, the magnetization m and the magnetic susceptibility 
Xm are determined by the g = ri„ = Fourier com- 
ponent of the order parameter field, and do not depend 
on the finite-frequency behavior of the critical modes. 
Therefore, dynamical scaling does not enter the tem- 
perature dependence of m and directly. In other 
words, their behavior is co mpletely determined by the 
equation of state, Eq. (3.6), whose coefficients acquire 
finite temperature corrections acco rdin g to Eqs. ( 2.10b| ) 
and ( ^.llb ). The tm term in Eq. (3.6) has a correction 
of the form m{p + TY''^ r^m{m + Tf-^. Similarly, the 
higher order terms t'np'^^^"" m"~^ have corrections pro- 



portional to {p + r)'^+i-" m" 



(m + T) 



d+l—n , 



This is the ballistic temperature scale discussed in Sec. 
IV A above. In addition, there is the fluctuation tem- 



perature scale that results fr om t he tempe ratur e depen- 
dence of t. By using Eq. ( 3.5a ) in Eq. ([4.1d| ) we see 
that [Tjbaii > [2^]fluct5 E^nd thus is the dominant temper- 
ature scale, for 1 < c? < 2, and for d > 5, but that for 
2 < d < 5, [Tjfluct is the relevant scale. All of these cor- 
rections can be summarized in the following scaling law 
for the magnetization m, 

m{t, T, H) = b-f^'^mith^'^ Th'l''\ Hb^^'") , (4.2a) 

where b is an arbitrary scale factor. The crossover expo- 
nent 



u = l/{d- 1) 
d/2{d~l) , 
3/(d+l) , 
J. = 1/2 , 



for 1 < < 2 
for 2 < d < 3 
for 3 < d < 5 
for d > 5 , 



(4.2b) 



describes the crossover from the quantum critical region 
to a regime whose behaviac, is, dominated by the clas- 
sical Gaussian fixed point. 0'E2l Note that, for the rea- 
sons explained above, the crossover exponent is not given 
by zv which one would expect from dynamical scaling. 
Also note the complicated behavior of the crossover ex- 
ponent as a function of the dimensionality, which is 
brought about by the competition between the ballistic 
and the fluctuation temperature scales. The result (j) = v 
that was reported in Refs. || and 20 was corre ct on ly for 
1 < d < 2 and d > 5. By differentiating Eq. ( [f.2a| ) with 
respect to H, we obtain the analogous homogeneity law 
for the magnetic susceptibility Xm, 

Xm{t, T, H) = f/'^x^ith'/^Tb^^'', Hb'^''') . (4.3) 



C. The specific lieat 

The scaling behavior of the specific heat is determined 
by the sum of the mean-field and the Gaussian fluctua- 
tion contribution to the free energy density /. The mean- 
field contribution follows immediately from Eq. (3.7). 



The Gaussian fluctuation contribution, /g, which gives 
the leading nonanalyticity for the specific heat at the 
critical point, can be calculated in conmlete analogy to 
the case of classical (j)^-i\ieoTy in d > 4.c^ Neglecting an 
uninteresting constant contribution to /g, we obtain. 



fa 



T 
2V 



E 



2ln{H/m + ad-iq'^^^ 
+02 + an \i^n\/q) 



\n[xdH/m- {xd - l)t + Od-iq 
+02 q^ + an \^n\/q) 



(4.4) 



Here Xd = d for 1 < d < 3, and = 3 for d > 3. The 
specific heat coefficient 7y is conventionally defined by. 



7v = cv/T = -d'f/dT' 



(4.5) 



Again we are interested only in scaling properties and 
not in exact coefficients. Keeping this in mind, an ade- 
quate representation of Eqs. (4.4) and (4.5) is given by 
the integral. 



7v 



dq 



T + q'i + q^ + Hq/r 



(4.6) 



with A an ultraviolet cutoff. 

Let us point out two interesting features of this result. 
First, for all dimensionalities in the range 1 < d < 3, 
"fv is logarithmically singular for T, TJ — > 0. This can be 
seen most easily fr om Eg. (|4.6|), and it is also true for the 
exact result, Eqs. (4^), ( [4.4| ). Such a d-independent loga- 
rithmic singularity is somewhat unusual. Wegner has dis- 
cussed how logarithmic corrections to scaling arise if |a-set 
of scale dimensions fulfills some resonance condition.c3 In 
the present case the relevant relation is that the scale di- 
mension of the free energy, d + z — 2z for 1 < d < 3, is 
a multiple of the scale dimension of T, which is z = d 
in this region. The fact that the logarithm appears in a 
range of dimensions, rather than only for a special value 
of d, is due to the dynamical exponent being exactly d in 
that range. Second, as discussed in Se c. |IV A| , two differ- 
ent temperature scales appear in Eq. (|4.(^). The first two 
terms in the denominator indicate that T ^ C as one 
would expect from dynamical scaling. However, the last 
term in Eq. (4.6) contains the magnetization, which in 
turn depends on [Tjbaii and [Tjfiuct- These two features 
imply that the scaling equation for jv should be written. 



-fvit,T,H) = e(3-d) In 6 

+F^{tb^l'',Tb^,Hb^'^l'',Tl 



(4.7) 



Note that the scale dimension of 7y , ignoring the loga- 
rithm, is zero in all dimensions. Since z > (jj/u for all 
d > 1, one can formally ignore the fourth entry in the 
scaling function since it is subleading compared to the 
second entry and its effects can be considered as 'correc- 
tions to scaling'. We emphasize that in contrast to the 
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magnetization and the magnetic susceptibility, the spe- 
cific heat does depend on the critical modes, and hence 
cont ains t he critical temperature scale. As mentioned in 
Sec. [V A , the latter is dominant when it is present, and 
provides an example for that. 



D. The disordered case revisited 

Let us finally rjficonsider the case of systems with 
quenched disorder.cl We do this partly to point out the 
remarkable analogy between the clean and dirty cases, 
and partly to correct a mistake in the results of Ref. ^ 
for dimensions d > 3. 

In the presence of disorder, the temperature or fre- 
quency far from criticality scales like the wavenumber 
squared. The baUistic temperature scale of Eq. ( 4.1a ) 
therefore gets replaced by a diffusive one, 



m 



diff 



(4.8) 



The s cale dim ensio ns for the other temperature scales, 
Eqs. ( pH ) - ( p^ ) remain valid, but the values of the 
exponents z, and 9 change. In the action, the disor- 
der leads to terms in add ition to, and structurally dif- 
ferent from, those in Eq. ( ^.7a ).a In particular, at order 
a second term appears whose coupling constant was 
denoted by U4 in Ref. ^. This operator, whose scale di- 
mension is [W4] — — — 4|, is the least irrelevant operator 
that is dangerous with respect to i, and one therefore 
has 9 = — \d— 4| in the fluctuation temperature scale, 
Eq. (4. Id). With z and v as determined in Ref. ^, this 



leads to a crossover exponent 

{2v = 2/{d-2) 
d/2 , 
4/(d-2) , 
2i^ = l , 



for 2 < d < VE+l 

for ^/5 + 1 < d < 4 

for 4 < rf < 6 , 

for d> 6 , 



(4.9) 



The result (j) — 2u of Refs. ^, was thus not correct for 
the (unphysical) dimensionality range V5+I < d < 6. 
The scaling behavior of the magnetiz ation and the mag- 
netic susceptibility is given by Eqs. (4.2a) and (4.2), re- 
spectively, with (f) from Eq. (|4.9|) , and all other exponents 
as given in Ref. ^ 



V. HERTZ'S MODEL REVISITED 

In this section we reexamine Hertz's original model of 
the ferromagnetic quantum phase transition. El We show 
that at tree level, the LGW theory for this model breaks 
down for related, but somewhat different reasons than in 
the realistic model above. Moreover, starting at one-loop 
order the renormalization group generates terms that are 
not in the bare action. As a result, the critical behavior 
of this model in d > 1 is not mean-field like, but rather 



the same as that of the more realistic model we have 
studied so far. 

Hertz's model differs in two respects fr om the more re- 
alistic one given by Eqs. (2.2) through (2.5). First, the 
interaction part of the action contains only the particle- 
hole spin-triplet channel that is decoupled in the deriva- 
tion of the LGW functional. Consequentely, the reference 
ensemble Sq consists of noninteracting electrons. Second, 
the coefficients x^"'' of the LGW functional are taken to 
be the correlation functions of a 3-d Fermi gas, irrespec- 
tive of the dimensionality of the space the magnetization 
is confined to. For this model the spin susceptibility Xs 
of the r efer ence system is simply the Lindhardt function, 
so Eq. (2.9) gets replaced by 



Xs(q, fin) = Co - C2 9^ - \ fln\/q 



(5.1) 



In comparison to the ana logo us ex pressio n for an inter- 
acting Fermi liquid, Eqs. (2^) and ( 2.10a ), the term pro- 
portional to g**"^ is missing, and Xs(q, = 0) is now 
analytic at q = 0. The resulting Gaussian part of the 
bare LGW functional has the form 



$2[M] ^\Y1 ^(^' + ^29' + an \nn\/q] 

xM(-q,-ri„) . (5.2) 



This action allows for a Gaussian fixed point with meaii= 
field static exponents and a dynamical exponent z = 3.EI 
Let us now investigate the stability of this Gaussian 
fixed point. At tree level, we can do this by calculat- 
ing the scale dimensions of the coefficients of the higher 
order terms. Hertz considered the higher order corre- 
lation functions x^"^ only in the limit il„ = 0, q — > 
where they are finite numbers. T he us ual power count- 
ing arguments in analogy to Sec. 
those terms are irrelevant for d 



[HB 



> 



show that all of 
1. The quartic one 
becomes marginal in d = 1 and relevant for d < I. This 
changes, however, if one considers the first order term in 
an expansion of the x'"'' for non-interacting electrons in 
powers of frequency. The reason is that, at nonzero ex- 
ternal frequency, frequency mixing effects occur that are 
similar to those brought about by an electron-electron 
interaction (which causes them even at zero external fre- 
quency). For power counting purposes, i.e. to determine 
the scale dimensions of the coefficients, it is not neces- 
sary to calculate the x*-"'' completely as functions of rt — 1 
frequencies and wavevectors. Guided by more complete 
calculations for the cases n = 3,4, 5, we have concluded 
that for power counting purposes it suffices to consider 
one independent frequency, and two directionally in- 
dependent wavevectors of equal length, |qi| = |q2| = |q|, 
that form an angle 9. With these simplifications the lead- 
ing term of the general coefficient is easily calculated. We 
find that for odd n the linear-in-frequency term in the ac- 
tion can symbolically be written 



Wn / dx W{x) n/q''- 



(5.3a) 
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with Wn some coupling constant, while the correspond- 
ing terms for even n are less relevant. The sa me power 
counting arguments that we used in Sec. Ill B| show that 
the scale dimension of Wn is [wn] = —{d~ l){n ~ 2)/2, 
which is negative for all d > 1. The coupling constants 
Wn arc thus irrelevant for all physical dimensions d > 2, 
and s eem to become marginal in d — 1. However, Eq. 
( ^.3a ) holds only if the two independent wavevectors are 
neither parallel nor antiparallel, i.e. it holds only if d > 1. 
For = 0, TT one finds a stronger singularity for the terms 
with even n, 



(a) S= ---O-- + 



■<i: 



w„ / dx Ar{x) n/q 



2n-3 



(5.3b) 



while the terms with odd n are less relevant. In dimen- 
sions d > 1 the parallel wavevectors form a set of measure 
zero and this stronger s ingul arity is of no consequence. 
In d = 1, however, Eq. (5.31) represents the generic be- 
havior of the terms of 0(M"). Power counting yields 
the scale dimensions to be [wn] — n — 2 in d — 1, so 
all of these coefficients are relevant operators. This is 
sufficient to conclude that the upper critical dimension 
is not one, but rather that the 1-d sytem is below its 
upper critical dimension, and will show nontrivial criti- 
cal behavior. This provides a technical explanation for 
Sachdev's observationB that Hertz's results in d < 1 can- 
not be correct. 

Moreover, the renormalization of the model beyond the 
tree level qu alita tively changes the form of the Gaussian 
action, Eq. (5^). Consider, for instance, the one- loop 
renormalization of <I>2 by the terms of order M'^ and M"*. 
The corresponding diagrams are shown in Fig. |l|. It is 
easy so see that these diagrams are equivalent to those 
that determine the spin susceptibility of intpcacting elec- 
trons, and have been calculated in Rcf. |l2|.E3 Renormal- 
ization therefore leads to a term proportional to q"^^^ in 
t he G aussian action, and hence to a $2 as given by Eq. 
la). We conclude that the critical behavior of Hertz's 



(1 



model for d > 1 is not mean-field like, but rather the same 
as that of the more realistic model dicussed in Sees. || - 

fvl. 



VI. DISCUSSION 

In this paper we have shown that clean itinerant quan- 
tum ferromagnets at zero temperature do not show, as 
was previously thought, uninteresting mean-field critical 
behavior for all dimensionalities d > 1. Rather, there are 
two upper critical dimensions. The first one is djj" = 1, 
above which the critical behavior is described by a Gaus- 
sian field theory but is not mean-field like, and the sec- 
ond one is dji""*" = 3, above which one does find mean-field 
critical behavior. The reason for this unusual behavior is 
soft modes that lead to an effective long-range interaction 
between the order parameter fluctuations. As is the case 
for classical models with long-range interactionspj this 



(b) 



■<:>■ 



n 



<>■■■■>■ 




FIG. 1. (a) Symbolic representation of the action. The 
dashed lines represent the external legs, and the polygons 
represent the vertices, (b) Renormalization of the two-point 
vertex to one-loop order. 



leads to nontrivial critical behavior that still can be de- 
termined exactly, djj""*" is the marginal dimension where 
the soft mode induced long-range interaction coincides 
with the usual r~('^+^) behavior. In this final section we 
discuss a few aspects of these results that have not been 
covered yet. 

First of all, both our approach and our results are 
remarkably similar to a recent treatment of disordered 
itinerant quantum ferromagnets.u In these papers it was 
found that the disorder induced diffusive excitations in a 
Fermi system with quenched disorder lead to similar, but 
stronger effects, with a long-ranged interaction between 
order parameter fluctuations that falls off like r~^^'*~^-', 
and three 'upper critical dimensionalities', viz. 2, 6, and 
4. The first two are analogous to d+ and d++ above. 
In addition, the critical exponents v, rj, and z take on 
their mean-field values for d > 4, while /3 and 5 remain 
anomalous between d = 4 and d — 6. In the present 
case, a different structure of the dangerous irrelevant 
variables makes the upper two special dimensions coin- 
cide. The same difference in the dangerous irrelevant 
variables leads to a difference in the temperature depen- 
dence of the equation of state as a function of dimension- 
ality: In the disordered case, the diffusive temperature 
scale with [rjdis = 2, which is analogous to the ballistic 
temperature scale [Tjbaii = 1 in the clean case, is domi- 
nant in the physical dimension d = 3. In the clean case, 
in contrast, the ballistic scale is subdominant compared 
to the fluctuation temperature scale for 2 < d < 5. 

The present paper shows that the basic concepts of 
Refs. ^ are not restricted to disordered systems. Indeed, 
an attempt to construct an effective field theory entirely 
in terms of the order parameter field for any phase transi- 
tion will break down (in the sense that it is impossible to 
construct a local effective theory) if there are soft or slow 
modes other than the order parameter fluctuations that 
couple to the order parameter. In the present case, the 
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spin-triplet particle-hole excitations that always exist in 
an interacting Fermi system, and that are distinct from 
the order parameter mode in that they are soft even in 
the paramagnetic phase, are such additional soft modes. 
They lead to the nonanalytic behavior of spi n de nsity 
correl ation functions that is displayed in Eqs. ( 2.10| ) and 
( ^.11 ), and hence to the effective long-ranged interaction 
between the order parameter modes. In general, the ap- 
pearance of such additional soft modes would call for the 
derivation of a different effective theory that does not in- 
tegrate out as many degrees of freedom, and that keeps 
all of the soft modes on equal footing. However, as in 
Ref. ^ we have opted here to work with the order param- 
eter effective theory after all, since it turns out that the 
difficulties introduced by the nonlocal character of the 
field theory can be overcome. Nevertheless it would be 
interesting to treat the same problem by means of a local 
theory that keeps more degrees of freedom explicitly. 

It should also be pointed out that our results depend 
crucially on the fact that the electronic spin density is 
a conserved quantity. If this feature was lost, e.g. due 
to some type of spin flip process, then the spin-triplet 
particle-hole excitations would acquire a mass or energy 
gap, and at scales larger than this mass the effective inter- 
actions between the order parameter fluctuations would 
be of short range. The asymptotic critical phenomena 
would then be described by a local order parameter field 
theory with mean-field critical behavior in all physical 
dimensions. At this point one might wonder whether 
the magnetization in the ordered phase, and magnetic 
fluctuations in the disordered one, act effectively as mag- 
netic impurities, and why this does not lead to an energy 
gap that invalidates our conclusions. The answer it that 
this effect has already been taken into account. In the 
ordered phase, the magnetization indeed acts as a cut- 
off, as has been discussed in connection with Eq. (3.6), 
and this leads to the nonanalyticity in the equation of 
state. In the disordered phase, the cutoff enters only via 
fluctuations, which are RG irrelevant with respect to the 
Gaussian flxed point. The effect therefore manifests it- 
self only in the corrections to scaling, not in the leading 
scaling behavior. 

We also mention that all of the qualitative points dis- 
cussed in Refs. ^ and ^ that had to do with the fact 
that one works above an upper critical dimension, ap- 
ply here as well. In particular, the presence of danger- 
ous irrelevant j/ariables means that some of Sachdev's 
general resultaJ are not directly applicable to the tran- 
sition discussed here. For instance, the Wilson ratio, 
W = {m/H)/{Cv/T), diverges at criticality, as it does in 
the disordered case, rather than being a universal num- 
ber as would be the case in the presence of hyperscaling. 
However, due to the different structure of the dangerous 
irrelevant variables that was already mentioned above in 
connection with the multiple upper critical dimensional- 
ities, some details are different between the disordered 
and clean cases. For instance, the scaling function F^, 
Eq. (4.7), is a function of T/H (if one neglects the sub- 



leading ballistic temperature scale) in agreement with the 
prediction of Ref. |[ while in the presence of disorder this 
is not the case.a 

We add one more remark concerning Hertz's original 
model that was discussed in Sec. ^ Our conclusion that 
a proper renormalization of that model leads to a crit- 
ical behavior that is the same as that of the realistic 
model solves the following paradox that would otherwise 
arise: Suppose one divided the interaction term in Hertz's 
model into two structurally identical pieces, one carrying 
some fraction of the interaction amplitude Ft, and the 
other the rest. Suppose one then applied the Hubbard- 
Stratonovich decoupling only to one of these pieces, and 
lumped the other into the reference ensemble. Then the 
reference ensemble would contain a nonzero Ft, and ac- 
cording to Ref. U2 the Gaussian action would contain 
the nonanalytic q^~^ term that leads to non-mean field 
critical behavior for all c? < 3. If one decouples all of 
the interaction term, on the other hand, then the ref- 
erence ensemble has Ft = 0, and in the absence of any 
other interaction amplitudes there is no way to generate 
a nonzero Fj by renormalization. If Hertz's model had 
indeed a critical behavior that is different from that of 
the realistic model, then the inevitable conclusion would 
be that the result depends on how exactly one performs 
the decoupling of the interaction term, which would be 
physically absurd. 

Finally, it should be mentioned that an experimen- 
tal corroboration or refutation of our results is proba- 
bly harder for the clean case discussed in the present 
paper than for the disordered case treated in Refs. |^. 
There are several reasons for this. First of all, the zero- 
temperature behavior can of course not be observed di- 
rectly. An experimental study would therefore have to 
concentrate on the crossover from the quantum critical 
behavior to the classical one that will occur if the classical 
transition point is at a sufficiently low temperature for 
the crossover point to be within the critical region. This 
requires a ferromagnet with as low a Curie temperature 
as possible. In addition, the parameter t that measures 
the distance from the critical point is, in the quantum 
case, not the temperature but rather the exchange inter- 
action or some other microscopic parameter that is hard 
to control. Both of these difficulties can be overcome rel- 
atively easily in the disordered case, where a change in 
the composition of a ferromagnetic alloy changes both 
the classical transition temperature and t. The quantum 
critical point is reached in the low-temperature limit in 
a sample whose composition is such that the Curie tem- 
perature just vanishes. As was discussed in Ref. U, this 
provides a manageable handle on t that has no obvious 
analog in the clean case. Furthermore, the differences 
between the quantum critical exponents in c? = 3 and the 
classical Heisenberg exponents are more pronounced in 
the disordered case than in the clean one. As we have 
seen, the 3-d critical behavior in the latter case is mean- 
field with logarithmic corrections to scaling. The loga- 
rithms would be hard to observe in any case, and the 
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p + q> J 



p + q, j k + q, j 



FIG. 3. Construction of a particle-hole spin-triplet vertex 
out of particle-particle vertices. 



p + qj 



k, i 



r3 = r. 



FIG. 2. The three basic interaction vertices. Straight lines 
with arrows denote particles and holes, wavy lines denote the 
interaction. The notation is the same as used in the text. 



mean-field exponents are much closer to 3-d Heisenberg 
exponents than those obtained in Ref. ^ 
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APPENDIX B: LOGARITHMIC CORRECTIONS 
TO SCALING FOR D = 3, AND FOR 1< D < 3 

There are three distinct mechanisms that produce log- 
arithmic corrections to scaling: (1) Marginal operators, 
(2) Wegner resonance conditions between a set of scale 
dimensions, and (3) logarithmic corrections to the scale 
dimension of a dangerous irreleicant operator. The first 
two mechanism are well known.Ej The third is operative 
only above an upper critical dimension, and is therefore 
of particular interest for quantum phase transitions. 

In the present case, logarithmic corrections to scaling 
arise due to mechanisms (2) and (3). Mechanism (2) 
produces corrections to the scaling of the specific h eat in 
all dimensions 1 < c? < 3, as was discussed in Sec. IV C . 
The third mechanism pr oduces corrections to scaling in 
d — 3. According to Eq. (2.11a), the correlation function 



X 



(4) 



contains a term \np in d = 3. Via Eq. ( p.q ) or 
( ^.7[ ) this leads, for instance, to a leading behavior of the 
spontaneous magnetization. 



m{t,H = 0) 



^1/2 



and at the critical point wc have 

771/3 



m{t = Q,H) 



Ml/H)] 



1/3 



(Bl) 



(B2) 



APPENDIX A: RENORMALIZATION OF THE 
PARTICLE-HOLE SPIN-TRIPLET 
INTERACTION VERTEX 



In this appendix we show that a nonvanishing particle- 
particle channel interaction, Fc 7^ 0, generates a particle- 
hole spin-triplet interaction Tt, even if Ff =0 in the bare 
system. 

Figure || shows the three basic types of interaction ver- 
tices. Now suppose that Ff — 0. Then a vertex of this 
type can be constructed by means of the diagrams shown 
in Fig. where all of the dashed lines represent Cooper 
channel interaction amplitudes. 
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